(N 
O 
O 
(N 



HUPD-0120, YITP-02-19, RCNP-Th02004 
hep-lat/0203025 

Heavy quark action on the anisotropic lattice 

Junpei Harada, 1 Hideo Matsufuru, 2 Tetsuya Onogi, 2 and Ayumu Sugita 3 

'Department of Physics, Hiroshima University, Higashi-hiroshima 739-8526, Japan 
2 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
'Research Center for Nuclear Physics, Osaka University, Ibaraki 567-0047, Japan 

(Dated: March 23, 2002) 

We investigate the O(a) improved quark action on anisotropic lattice as a potential framework for 
the heavy quark, which may enable precision computation of hadronic matrix elements of heavy-light 
mesons. The relativity relations of heavy-light mesons as well as of heavy quarkonium are examined 
on a quenched lattice with spatial lattice cutoff a& 1 ~ 1.6 GeV and the anisotropy £ = 4. We find 
that the bare anisotropy parameter tuned for the massless quark describes both the heavy-heavy 
and heavy-light mesons within 2% accuracy for the quark mass a a mQ < 0.8, which covers the charm 
quark mass. This bare anisotropy parameter also successfully describes the heavy-light mesons in 
the quark mass region a a m,Q < 1.2 within the same accuracy. Beyond this region, the discretization 
effects seem to grow gradually. The anisotropic lattice is expected to extend by a factor £ the quark 
mass region in which the parameters in the action tuned for the massless limit are applicable for 
heavy-light systems with well controlled systematic errors. 

PACS numbers: 12.38.Gc 
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I. INTRODUCTION 

Recent experimental progress in flavor physics to look 
for the effect of new physics requires precise theoretical 
predictions from the standard model. However, model 
independent calculation of hadronic matrix elements is 
difficult because of nonperturbative nature of QCD. The 
lattice QCD simulation is one of most promising ap- 
proaches in which the systematic uncertainties can be 
reduced systematically jl). The ultimate goal of this 
work is to construct a framework for lattice calculations 
of hadronic matrix elements in a few percent system- 
atic accuracy, as required by the experiments in progress. 
More practically, the level of accuracy requested from the 
CLEO-c experiment Q as well as from the B Factories 
H [| are about 2 %. This paper intends to propose a 
systematic program to achieve this accuracy making use 
of the anisotropic lattice QCD. 

In a lattice calculation of heavy quark systems such 
as the charm and bottom, one needs to control the large 
discretization error of 0((am,Q) n ). Extensive studies in 
various approaches have achieved much progress in un- 
derstanding the heavy quark systems. However, it still 
seems difficult to achieve aforementioned systematic ac- 
curacy with the current techniques for several reasons. 
Here we categorize them into three types and summarize 
their advantages and disadvantages from the viewpoint 
of precision study of matrix elements in flavor physics. 

1 ) Effective theories. One approach is to use descrip- 
tions of the heavy quark based on the heavy quark effec- 
tive theory using nonrelativistic QCD (H). As an advan- 
tage, the action incorporates higher order correction in 
spatial momenta and can easily remove the all the mass 
dependent errors at the tree or one-loop level. However, 
since the theory does not have the continuum limit, one 
cannot remove the discretization errors by extrapolating 
the results on lattices with finite spacings. Another dis- 



advantage is that the nonperturbative renormalization is 
difficult due to strong mass dependence. Precise mea- 
surements of the weak matrix elements for the heavy- 
light mesons below 10% are therefore difficult. 

2) Relativistic framework. The most straightforward 
approach is to use the 0(a) improved Wilson action in 
natural way for relatively lighter heavy quark, and to ex- 
trapolate the results in 1 jm according to the heavy quark 
effective theory. Since the theory has the continuum 
limit, the discretization errors can be removed by an ex- 
trapolation. The perturbative errors can also be avoided 
by employing the nonperturbative renormalization. In 
practice, however, lattice artifacts of 0((arriQ) 2 ) give 
dominant errors Q so that the precise calculation will 
remain difficult for the next few years even in quenched 
approximation. Brute force improvement by decreasing 
lattice spacing and employing large lattices quickly in- 
crease the simulation cost, and therefore not a realistic 
solution for the calculation of matrix elements. 

3) Fermilab approach. The Fermilab approach links 
the above two approaches @, || ■ In the heavy quark mass 
region, the 0(a) improved Wilson action with asymmet- 
ric parameters is reinterpreted as an effective theoreti- 
cal description just like the nonrelativistic QCD action. 
Since the action reduces to the conventional O(a) im- 
proved Wilson action for small masses arriQ <C 1, it has 
in principle smooth continuum limit. The disadvantage 
is that it is not known how to extrapolate the results on 
lattices with amg > 1, which are currently unavoidable 
in particular for systems with 6-quark. For such a quark 
mass region, a mass dependent tuning of parameters in 
the action is required for proper improvement, although 
a systematic tuning beyond the perturbation theory is 
still a challenge Therefore precise calculation below 
10% level are nontrivial for the weak matrix elements of 
the heavy-light mesons. 

Therefore, it is desirable to develop a new framework 



2 



for the heavy quark with the following features, (i) The 
continuum limit can be taken; (ii) A systematic improve- 
ment program, such as the nonperturbative renormaliza- 
tion technique [fiof , can be applied not only to the pa- 
rameters in the action but also to the operators; (hi) A 
modest size of computational cost is required for system- 
atic computation of matrix elements. 

The anisotropic lattice, on which the temporal lattice 
spacing a T is finer than the spatial one a a , is a candi- 
date of such a framework jn], [l^] . The approach is fun- 
damentally along the line of Fermilab approach on the 
anisotropic lattice. However, large temporal lattice cut- 
off is expected to improve the above problems drastically. 
The most crucial one is that the mass dependence of the 
parameters in the action may become so mild in the re- 
gion of practical interest that one can adopt the 0(a) 
improving clover coefficient determined in the nonper- 
turbative renormalization technique. On the other hand, 
standard size of spatial lattice spacing keeps a total com- 
putational cost modest. Therefore, extrapolation of the 
result to the continuum limit may be possible with keep- 
ing systematic uncertainties under control with sufficient 
accuracy. Whether these observations practically hold 
should be examined numerically, as well as in the pertur- 
bation theory. 

Our form of quark action on the anisotropic lattice is 
founded in Refs. [jjjl [y|, where the perturbative results 
for light and heavy quarks and the simulation results of 
the light quark were studied. In this paper, we focus on 
the heavy quark and study the mass dependence of the 
the Lorentz symmetry breaking effect in order to under- 
stand the mass range for which a consistent description 
is possible. This signals the applicability of anisotropic 
lattice to the heavy quark systems. We investigate the 
heavy-heavy and heavy-light mesons on a quenched lat- 
tice of anisotropy £ = 4 and spatial cutoff a" 1 ~ 1.6 
GeV. The heavy quark mass is varied from the charm 
quark mass to about 6 GeV to examine an applicability 
of the action to these systems. 

This paper is organized as follows. In the next sec- 
tion, we first give our quark action, which is discussed 
in detail in Ref. We then give our conjecture on 

the lattice spacing dependence of the anisotropy param- 
eter toward the continuum limit, and discuss the advan- 
tage of the anisotropic lattice compared to the isotropic 
one. In Section III, we observe the tree-level expecta- 



tion of the mass dependence of the 0(a 2 ) terms in the 
quark dispersion relation and study how the anisotropic 
parameters or the breaking effect of relativity behave as a 
functions of the heavy quark mass. Section IV describes 
the numerical simulation. We compute the heavy-light 
and heavy-heavy meson spectra and dispersion relations 
with two sets of heavy quark parameters. In one set (Set- 
I), the bare anisotropy is set to the value for the massless 
quark. The other set (Set-IT) adopt the result of the mass 
dependent tuning using the heavy-heavy meson disper- 
sion relation as obtained in Ref. |L3]]. We observe the 
mass dependence of the renormalized anisotropy in or- 



der to probe the breaking of relativity for heavy quark. 
We also observe how the inconsistency among the bind- 
ing energies of heavy-heavy, heavy-light and light-light 
mesons JlJ, ^H) grows as a function of mass. The last 
part of this section discusses the hyperfine splitting of 
the heavy-light meson. In Section |V| we summarize the 
result of simulation, and draw our perspective on further 
development of this framework, as our conclusion in this 
paper. 



II. FORMULATION 

A. Quark action 

We adopt the following quark action in the hopping 
parameter form Jl2|, : 



S F = ^2 i>{x)K{x, y)tp(y), 



(1) 



x.y 



Sf = S x . y -K T (1 - ^)Ui{x)8 x+ i y 

+ (1+74)^(^-4)^ 

i 

-K a c E cruFu{x)8 x>v - rn a c B a lj F lj (x)S x ^,(2) 

i i>j 

where k ct and the spatial and temporal hopping 

parameters, and are related to the bare quark mass and 
bare anisotropy as given below. The parameter r is the 
spatial Wilson coefficient, and the parameters ce and 
cb are the clover coefficients for the O(a) improvement. 
Although the explicit Lorentz symmetry is not manifest 
due to the anisotropy in lattice units, it can be restored in 
principle for physical observables in physical units at long 
distances up to errors of 0(a 2 ) by properly tuning k ct /k t , 
r, ce and cb for a given n a . The action is constructed 
in accord with the Fermilab approach and hence ap- 
plicable to an arbitrary quark mass, although a mass de- 
pendent tuning of parameters is difficult beyond the per- 
turbation theory. This may be circumvented by taking 
a^T 1 ^> niQ, with which the mass dependence of param- 
eters are expected to be small so that 0(a) Symanzik 
improvement program for the heavy quark can be ap- 
plied. To study whether this expectation really holds or 
not is the main subject of this paper. 

In present study, we vary only two parameters n a and 
k t with fixed other parameters. We set the Wilson pa- 
rameter as r = l/£ and the clover coefficients as the 
tadpole-improved tree-level values, ce = 1/ugU 2 ., and 
cb = The tadpole improvement fl7| is achieved 

by rescaling the link variable as Ui(x) — > Ui(x)/u a and 
Ua(x) — > Ua{x)/u t , with the mean-field values of the spa- 
tial and temporal link variables, u a and u T , respectively. 
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Instead of n a and k t , we introduce k and 7f(= 1/C) as 

11, , 

- = 2(7f + 3r - 4) = 2(m j F a + 4), 

7f = 1/C = K T u T /n a u a . (3) 

The former controls the bare quark mass and the latter 
corresponds to the bare anisotropy. 



B. Mass dependence of anisotropy parameter 

On an anisotropic lattice, one must tune the param- 
eters so that the hadronic states satisfy the relativistic 
dispersion relations. In general, the lattice dispersion re- 
lation for arbitrary values of £ is described in lattice units 
as 

E 2 {p) = m 2 +p 2 /e F + 0{p i ). (4) 

In the above expression the energy E and the rest mass m 
are in temporal lattice units while the momentum p is in 
spatial units. The parameter £f introduced in this equa- 
tion characterizes the anisotropy of the quark fields. The 
difference between the quark field anisotropy £f and the 
gauge field anisotropy £ probes the breaking of relativity. 
Therefore, the calibration is nothing but the tuning of 
£(= 1/7f) f° r a given k so that £f equals £ and hence 
the relativistic dispersion relations are satisfied. Let us 
call the tuned parameter as C*(= 1/7f)- C* depends on 
the quark mass and can in general be different from unity 
even for the case of isotropic lattice. The requirement of 
the relativity relation automatically enforces that the rest 
and kinetic masses are equal to each other. In this sense, 
our calibration procedure of anisotropic lattice action is 
a natural generalization of the Fermilab approach. 

Now we consider the tuning of anisotropic parameter 
either on isotropic or anisotropic lattice, for fixed physical 
quark mass. The anisotropy parameter £ = 1/t> can be 
tuned using the dispersion relation of either heavy-heavy 
or heavy-light meson in mass dependent way. Alterna- 
tively, one can also adopt the value tuned for the light- 
light mesons neglecting the mass dependence. These 
procedures of calibration can give different results due 
to the discretization errors. We would like to know: 1) 
the value of lattice spacing a!^ h ~ hl above which the cali- 
brated parameters £* using the heavy-heavy and heavy- 
light mesons differ from each other by more than a cer- 
tain accuracy e acc (say, 2%), and 2) the value of lattice 
spacing a^ l ~ n above which the calibrated parameters £* 
using the heavy-light and light-light mesons are different 
by more than a certain accuracy e' acc . 

Naively speaking, the difference between £*'s with 
heavy-heavy and heavy-light mesons originates from the 
0((a a p) 2 ) effects, where p is the typical quark momenta 
in heavy-heavy meson, p ~ urtiQ. On the other hand, 
the deviation of £* with heavy-light mesons from the 
£* with light-light mesons is due to the 0((a T rriQ) 2 ) = 



0((a a mQ/^) 2 ) effects. With the criteria that these errors 
stay within the required accuracies, we obtain 

a h a h - hl ~ V^c/(am Q ), 
4 l ~ U ~ V^TcH/mQ. (5) 

Figure |l| conjectures the expected behavior of the 
tuned anisotropy parameter £* for the isotropic (£ = 1) 
and anisotropic (£ = 4) lattices. The tuned anisotropy 
normalized by the value at massless limit is displayed for 
fixed physical quark mass, around charm quark mass, as 
an example. In the case of isotropic lattice, as shown 
in the left panel of Fig. Q, a^ l ~ u < a^ h ~ hl is expected 
to hold. In this case, for a sufficiently small lattice 
spacing a a < a£ , both the heavy- heavy and heavy- 
light mesons can be successfully described by setting the 
anisotropy parameter to the value in the massless limit. 
For a coarser lattice spacing a^ l ~ u < a a < a'^ h ~ hl , mass 
dependent tuning, namely the Fermilab approach, is nec- 
essary, but the description of both the heavy-heavy and 
heavy-light mesons still works. For an even coarser lattice 
spacing a'^ h ~ hl < a a , simultaneously consistent descrip- 
tion of heavy-heavy and heavy-light mesons no longer 
works, because of severe discretization effect in the heavy 
quarkonia. 

Now let us consider the anisotropic lattice case with 
£ = 4. The expected mass dependence of £* is schemat- 
ically represented in the right panel of Fig. |l|. As is ob- 
vious in Eq. (0) and will be studied in the following sec- 
tions, the anisotropy does not improve the situation for 
heavy-heavy meson. Therefore a^ h ~ hl remains roughly 
of the same size as in the isotropic lattice case. On the 
other hand, Eq. (|J) and the studies in the following sec- 
tions show that for £ = 4, a^ l ~ u is expected to increase 
and a^ h ~ hl < a^ l ~ u occasionally holds. This latter sit- 
uation is particularly probable for quark with not very 
large mass, and expected to hold up to the bottom quark 
mass region. In this case, for sufficiently small lattice 
spacing a a < a'^ h ~ hl , both the heavy-heavy and heavy- 
light mesons can be successfully described by setting the 
anisotropy parameter to the value in the massless limit. 
For a coarser lattice spacing a^ h ~ hl < a a < a^ l ~ u consis- 
tent description of heavy-heavy and heavy-light mesons 
no longer works. Nevertheless the heavy- light meson 
is successfully described with the anisotropy parameter 
tuned for the massless limit. Finally for an even coarser 
lattice spacing < a a , the genuine Fermilab ap- 

proach, namely mass dependent tuning of £, is required 
to describe the heavy-light mesons. 

In the range of lattice spacing a a < a™' , if we give up 
the description of heavy quarkonia and only consider the 
heavy-light systems, one can successfully adopt the tuned 
value of £ at the massless limit. This is the most impor- 
tant region for our approach, since there the mass depen- 
dence of parameters in the action other than £ = 1/7f, 
such as cf and eg, are also expected to be small, and 
hence the result of tuning at the massless limit is also 
considered to be valid. For the tuning of clover coefficient 
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1 1 

a, [GeV 1 ] a a [GeV 1 ] 

FIG. 1: The conjecture of lattice spacing dependence of the anisotropy parameter for fixed quark mass, for isotropic (left) 
and anisotropic (right) lattices. Horizontal lines with arrows roughly represent the region a a < a„ , while vertical thick lines 
correspond to a^ h ~ h , above which the f*'s determined from heavy-heavy and heavy-light mesons no longer equal. In this figure, 
a heavy quark with roughly the charm quark mass is considered as an example. 



for the massless quark, it may be possible to apply the 
nonperturbative renormalization technique [ fl0| , which is 
currently most efficient procedure to remove the 0(a) ef- 
fect beyond perturbation theory. The most advantage of 
the anisotropic lattice is that this region a < a% l ~ u is ex- 
pected to be extended by a factor £ as is found in Eq. (||), 
so as to cover the heavy quark mass region of practical 
interests. Therefore the anisotropic lattice is a candidate 
of framework for heavy-light systems, which satisfies the 
required characters (i), (ii), and (iii) mentioned in the 
Introduction. Also for heavy-heavy systems, it would be 
important to study quantitatively the heavy quark mass 
at which the action starts to fail in correct description. 

This paper is dedicated to a numerical study whether 
these expectations actually holds. Using Eq. (g), similar 
expectation can be deduced for the heavy quark mass de- 
pendence with a fixed lattice spacing. In this paper, since 
we work with only one lattice spacing, instead of study- 
ing the lattice spacing dependence with a fixed heavy 
quark mass, we study the mass dependence with a fixed 
lattice spacing. This means to specify in which region 
considered above the present a a is located for each value 
of quark mass. This subject is discussed in more detail 
with the tree level analysis in the next section, and with 
the numerical simulation in the successive section. 



III. EXPECTATION FROM TREE LEVEL 
ANALYSIS 



In this section, we review the tree level analysis of 
the heavy quark action in the Fermilab formulation. In 
the following, we express the energy and the momenta 
in physical units so that the dependence on the lattice 
spacings a T and a a is explicitly shown. The dispersion 



relation of the heavy quark on anisotropic lattice is given 
as % H 



cosh(a T i?) = 1 + 

[Mp + 2<E t sm 2 {a aPl /2)} 2 + C 2 sm 2 (a aPi ) 
2[1 + m + 2r( sm 2 {a aPl /2)} 



(6) 



where mo = a T mo . This leads to the following dispersion 
relation. 



+A 1 al(p'f + A 2 a',Y,Pt 

i 

where M\ is the pole mass, 

Mi = log(l + m )a 7 r 1 > 
and the anisotropy of the quark field at the tree level 



(7) 



(8) 



= log(l + m )£ 2 



2C 2 < 



mo(2 + mo) 1 + m 



The third and fourth terms in Eq (|7|) represent the 
0((ap) 2 ) errors, where Ai and A 2 are given as 



£ 2 



2C 2 



rC 



m (2 + TO ) 1 + mo 
8C 4 



log(l + m ) , 

Vmg(2 + TO ) 3 

4C 3 [C + 2r(l+m )] r 2 ( 2 



77*2(2 + m ) 2 (1+mo) 2 



(10) 
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FIG. 2: The mass dependence of tuned at the tree level 
for the isotropic lattice and anisotropic lattice with £ = 4. 
The horizontal axis is the kinetic quark mass in spatial lattice 
units. 



0.0 1.0 2.0 3.0 4.0 

FIG. 3: The kinetic mass dependence of £f/£ at the tree level. 
We do not show the the result of mass dependent tuning of £ 
since £f/£ = 1 is satisfied by definition. 



M = -ylog(l 



mo) 



2C 



m {2 + m ) 4(1 + mo) 

(11) 

The above coefficients are derived using the Mi, M%, M4 
and u>4 in Ref. extended to anisotropic lattice by re- 
placing certain parameters as given in Ref. [ jl2| . 

Fig. H shows the mass dependence of f tuned by re- 
quiring £p = £, both on the isotropic lattice and on the 
anisotropic lattice with £ = 4. In the latter case, it is 
clear that the mass dependence is drastically reduced so 
that taking the value of £ with massless limit is a good 
approximation over a wide range of quark mass, in con- 
trast to the case of the isotropic lattice. 

In the following, let us consider either the case with 
mass dependent tuning of £ (denoted by "Fermilab" in 
the figures) or with the £ at the massless limit ("mass- 
less" ) . Figure || shows the breaking effect of relativity re- 
lation, £f/£, according to Eq. (^). The kinetic mass M 2 
of quark is related to £f/£ as M 2 = A/i£f7£ 2 - F° r the 
mass dependently tuned C, £f trivially equals £. Without 
tuning, the mass dependence of £p is drastically reduced 
in the case of £ = 4, as is expected from that the de- 
viation of £f/£ from unity is the 0{a T m,Q) error. In 
Figure ^, one find that the coefficient A\ , whose devia- 
tion from zero signals the 0(o?) effect, is also drastically 
reduced on the anisotropic lattice, either with or without 
mass dependent tuning of £. On the other hand, 0(a 2 ) 
error from A 2 is not reduced but somewhat worse on the 
anisotropic lattice as is seen in Figure |j| However, the 
coefficient A 2 is not larger than the value at the mass- 
less limit. As a general feature for all of the £f/£, A\ 
and A 2 , the mass dependences are drastically reduced 
on the anisotropic lattice. On the other hand, mass de- 
pendent tuning of £ on isotropic lattice does not reduce 
the coefficients A\ and A 2 , while it completely removes 
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q=l, massless 

|=4, massless 

£=1, Fermilab 

^=4, Fermilab 



0.0 



1.0 



2.0 

a„ M, 



3.0 



4.0 



FIG. 4: The kinetic mass dependence of the 0((ap) 2 ) coef- 
ficient Ai in the dispersion relation (Q). "massless " denotes 
the result of tuning of £ in the massless limit. "Fermilab" 
denotes the result of mass dependent tuning of £ by requiring 

£f = C 



the discrepancy between £p and £. In order to reduce 
these errors one has to introduce higher order terms as 
is pointed out in Ref. (7) . 

From these results, one can expect the following about 
the error of the heavy quark systems. For the heavy- 
light systems, since heavy quark momenta typically takes 
the size of the hadronic scale, Aqcd — 200-500 MeV, 
the 0((ap) 2 ) errors are well under control even for large 
a a M 2 , if one keep a^ 1 ^> Aqcd- While the largest error 
originates from A 2 , this is of the same order as that in 
the light quark systems. For the heavy- heavy system, 
the situation is quite different. Since the heavy quark 
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FIG. 5: The kinetic mass dependence of the 0((ap) 2 ) coef- 
ficient A2 in the dispersion relation (pi). 



momenta are typically p ~ arriQ, The 0{{ap) 2 ) errors are 
expected to be as large as Ai(a.a a mQ) 2 , A2(aa a mo) 2 . 
Again the A^ gives the largest contribution, but the size 
of the error is uncontrolled when a a m,Q is large. 

To summarize, the anisotropic lattice largely reduces 
the discretization effects represented by £f/£~ 1 and Ai, 
while it does not improve the Ai . For the heavy- light sys- 
tems, this suffices for a computation with discretization 
effects under control. On the other hand, when a a M 2 is 
of order of unity, the anisotropic lattice does not improve 
the situation for heavy quarkonia, because of severe ef- 
fect of A2 in these systems. Although mass dependent 
tuning of Q further removes the deviation of £p from £, 
the C tuned for massless quark also stay a good approxi- 
mation if a T mo <C 1 holds. This is the most advantage of 
anisotropic lattice compared with the isotropic Fermilab 
approach, in which the mass dependence of the param- 
eter is much stronger. These observations are in accord 
with our conjecture in the last section for developing the 
anisotropic Fermilab formulation. 



IV. NUMERICAL SIMULATION 

This section numerically examines the idea described 
in previous sections. For this purpose, we perform sim- 
ulations with two series of heavy quark parameters. For 
Set-I the anisotropy parameter is set to the value at the 
chiral limit, jpirriq = 0), just as same as for the light 
quark. The Set-II adopts the fully tuned anisotropy pa- 
rameter 7^ using heavy-heavy mesons. This calibration 
is done in the first subsection. For the heavy-heavy and 
heavy-light mesons, the rest and kinetic masses are ob- 
tained with two sets of parameters. Two quantities are 
used to probe the breaking effect of the relativity: the 
fermionic anisotropy £p for heavy-heavy and heavy-light 



FIG. 6: The result of calibration in the heavy quark mass 
region (filled circles), together with the previous result in 
Ref. 13 . The fits are performed with all available data in- 
cluding previous results. 

mesons, and the inconsistency among the binding ener- 
gies of heavy-heavy, heavy-light, and light-light mesons. 
These behaviors give us a hint to judge the regions in 
which our framework can be applied. The last subsec- 
tion treats the heavy-light meson spectrum, and observe 
how the hyperfine splitting behaves with meson kinetic 
mass. 



A. Calibration in heavy quark region 

The simulation parameters used in this paper is the 
second set in Ref. the quenched anisotropic lattice 
of the size 16 3 x 128 generated with the standard pla- 
quette gauge action with (/?, 7g) = (5.95,3.1586), which 
correspond to the renormalized anisotropy £ = 4 [jlSfl and 
spatial lattice cutoff a~ 1 = 1.623(9) GeV set by hadronic 
radius rn |L9| . The mean- field values in the quark action 
are set to the mean values of link variables in the Landau 
gauge: u„ = 0.7917 and u T = 0.9891. 

In Ref. [ |i"3"| , the optimum bare anisotropy 7^ is de- 
termined using the dispersion relation of mesons with 
degenerate quark masses, and the resultant values of 7^ 
are well represented by a linear function 

4= co+c2 < m «=2i(]-£)' (12) 

where for the present lattice Co = 0.2490(8), (2 = 
0.189(15), and k c = 0.12592(6). 

In this paper, we use seven values of k for heavy quark 
(k/j), which roughly cover 1-6 GeV. Three of them al- 
ready appeared in Ref. jl3j . We start with the calibration 
of remaining four values of k in the heavy quark region in 
the same manner as in Ref. [fl3f . The values of k used are 
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TABLE I: Quark parameters and the result of calibration. The result for k — 0.1100, 0.1020, and 0.0930 are taken from 
Ref. The number of configurations is 200 for each value of the hopping parameter, except for k = 0.1100 for which the 

simulation is carried out with 300 configurations. 





m 


input 7i? 


*(PS) 


*(V) 


If 


5~fF 


0.1100 


0.1437 


3.90, 4.00 


3.945(28) 


3.946(36) 


3.946(29) 


-0.001(20) 


0.1020 


0.2328 


3.90, 4.00 


3.848(26) 


3.847(34) 


3.847(28) 


0.001(15) 


0.0930 


0.3514 


3.70, 3.80 


3.688(23) 


3.687(29) 


3.688(25) 


0.000(11) 


0.0840 


0.4954 


3.39, 3.48 


3.471(20) 


3.462(23) 


3.467(21) 


0.0092(60) 


0.0760 


0.6521 


3.04, 3.15 


3.205(20) 


3.191(23) 


3.199(21) 


0.0139(45) 


0.0700 


0.7930 


2.85, 2.95 


2.946(23) 


2.932(23) 


2.939(23) 


0.0134(52) 


0.0630 


0.9914 


2.50, 2,60 


2.584(22) 


2.561(24) 


2.573(23) 


0.0237(35) 



listed in Table | together with the result of calibration. 
The second column is the naive estimate of bare quark 
mass according to Eq. (|l2]). For the heaviest case, too 
is almost unity in temporal lattice units, and therefore 
the breaking effect of relativity may be visible. Here we 
note that for heavier quark masses the difference of 7^ 
for pseudoscalar and vector mesons, 



til 



*<V) *(PS) 

7f ~7_f , 



(13) 



is sizable beyond the statistical fluctuations. This signals 
that the quarkonium system is not properly described 
within the present framework at this lattice spacing. This 
problem will again be discussed later in terms of the 
fcrmionic anisotropics for heavy-heavy and heavy-light 
mesons. 

Figure || shows the result of calibration. The result 
is well fitted to a linear form in to^ or a quadratic form 
in m q , in spite of large quark mass. The fits including 
previous data in Ref. [jl3| result in 

Linear: Co = 0.2510(6), C2 = 0.1437(26). (14) 
Quadratic: Co = 0.2473(12), Ci = 0.0267(72), 

C2 = 0.1151(81). (15) 

In both cases, the values of Co is close to the value at 
the mean-field tree level, = 0.25. Since there is no 
reason that the small quark mass dependence persist up 
to m q ~ a~ , the difference between C2's in linear fits 
with present data and with only previous data in lighter 
quark mass region is irrelevant, and above result just 
shows that the quark mass dependence of 1/7^ is still 
tractable at this mass region. 



B. Relativity breaking effect 

In this subsection we compute the heavy-heavy and 
heavy-light meson spectra and the dispersion relations 
using two sets of heavy quark parameters. For the first 
one, Set-I, the bare anisotropy is set to the value of 7^ in 
the chiral limit obtained in Ref. namely jf = 4.016 
at f3 = 5.95, for all quark masses. The second one, Set-II, 



adopts the above result of mass dependent calibration. 
We use the same hopping parameter values Kh for the 
heavy quark as given in the previous subsection. For the 
light quark, we use single value, k; = 0.1235. The value 
of 7i? at this ki is set to the value in the chiral limit, 
as same as in the light hadron spectroscopy in Ref. Jl3| . 
We regard that at n = m the quark mass is sufficiently 
light for present purposes, and do not extrapolate to the 
chiral limit. The numbers of configurations used are 200 
and 500 for heavy- heavy and heavy- light meson masses, 
respectively. 



The meson spectra are listed in Table [II| and III for 
heavy- heavy and heavy-light mesons, respectively. In 
these tables, we also list the results of £f for each me- 
son channel, and difference between them, S^p = — 

^ < F PS \ If the anisotropic lattice action does not describe 
the quarks inside mesons in respecting the relativity re- 
lation, the breaking effect appears in the dispersion re- 
lations of mesons. Therefore the deviation of fermionic 
anisotropy £p from £ signals the breaking effect of rela- 
tivity. In the following, we first discuss on the result of 
£f for Set-I parameters, namely with 7^ tuned for mass- 
less limit, and then briefly summarize the result for Set-II 
parameters. 

Figure ^ displays the heavy quark mass dependence 
of £f for Set-I. The horizontal axis is the bare quark 
mass to ij in temporal lattice units. The behaviors of £p's 
are well in accord with the expectation in Sec. II. For 
quantitative discussion, let us consider the case that the 
required accuracies to define the a^. h ~ hl and a„ are 2 
%, namely e acc = e ac .c = 0.02. £p's from the heavy-heavy 
and heavy-light mesons disagree beyond this accuracy 
at TTiq > 0.2 (oo-toq > 0.8). Therefore one must keep 
vriq < 0.2 to avoid large systematic uncertainty in the 
heavy quarkonia. On the other hand, the £p from the 



heavy- light mesons are rather close to £, and 7^(7 



0) 



can be applied up to m q ~ 0.3 within presently required 
accuracy. For quark mass larger than this value, the dis- 
crepancy between £p 's from the pseudoscalar and vector 
mesons gradually grows beyond the statistical error. This 
signals the growth of systematic error. In the region of 
Trig < 0.3 such effect is sufficiently small. 
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TABLE II: Heavy-heavy meson spectrum for Set-I and Set-II parameters obtained with 200 configurations. 







m PS 


my 


my — mps 


SF 


£\ V ) 

SF 


A V I (\PO) 

SF Si? 


Set-I 0.1100 


4.016 


0.42468(23) 


0.43755(33) 


0.01288(18) 


4.069(34) 


4.068(46) 


-0.001(24) 


0.1020 


4.016 


0.58409(25) 


0.59358(34) 


0.00949(13) 


4.148(27) 


4.143(35) 


-0.005(16) 


0.0930 


4.016 


0.76947(24) 


0.77692(31) 


0.00745(10) 


4.292(22) 


4.294(28) 


0.002(10) 


0.0840 


4.016 


0.96746(25) 


0.97358(31) 


0.00612( 8) 


4.527(24) 


4.537(29) 


0.010( 8) 


0.0760 


4.016 


1.15894(27) 


1.16411(33) 


0.00517( 8) 


4.821(30) 


4.839(35) 


0.018(10) 


0.0700 


4.016 


1.31453(27) 


1.31896(31) 


0.00443( 7) 


5.146(33) 


5.187(38) 


0.041( 8) 


0.0630 


4.016 


1.51262(27) 


1.51617(31) 


0.00355( 5) 


5.650(39) 


5.726(43) 


0.076( 7) 


Set-II 0.1100 


3.946 


0.42942(23) 


0.44248(33) 


0.01306(18) 


4.005(29) 


4.003(41) 


-0.002(22) 


0.1020 


3.847 


0.60172(25) 


0.61146(33) 


0.00975(14) 


4.002(24) 


4.000(32) 


-0.003(14) 


0.0930 


3.688 


0.81797(24) 


0.82574(31) 


0.00777(10) 


3.995(19) 


4.000(24) 


0.005( 9) 


0.0840 


3.467 


1.07512(25) 


1.08167(30) 


0.00655( 7) 


3.996(18) 


4.005(22) 


0.009( 6) 


0.0760 


3.199 


1.35939(26) 


1.36509(31) 


0.00571( 7) 


4.002(22) 


4.011(26) 


0.010( 6) 


0.0700 


2.939 


1.62650(27) 


1.63164(31) 


0.00514( 6) 


3.993(23) 


4.009(26) 


0.016( 5) 


0.0630 


2.573 


2.02101(28) 


2.02550(31) 


0.00449( 5) 


3.988(23) 


4.013(26) 


0.025( 4) 



TABLE III: Heavy-light meson spectrum for Set-I and Set-II parameters obtained with 500 configurations. 







m PS 


my 


my — mps 


APS) 

SF 


Av) 

Sf 


Av) APS) 

SF Sf 


Set-I 


0.1100 


0.29108(27) 


0.30997(51) 


0.01889(39) 


3.985(36) 


4.006(49) 


0.020(44) 




0.1020 


0.37705(30) 


0.39106(50) 


0.01401(34) 


3.994(39) 


4.031(51) 


0.036(42) 




0.0930 


0.47595(36) 


0.48643(57) 


0.01048(36) 


4.032(51) 


4.101(73) 


0.070(55) 




0.0840 


0.58123(46) 


0.58123(46) 


0.00805(43) 


4.088(73) 


4.21(11) 


0.119(77) 




0.0760 


0.68279(55) 


0.68916(81) 


0.00638(44) 


4.148(99) 


4.30(15) 


0.147(87) 




0.0700 


0.76546(64) 


0.77079(89) 


0.00534(45) 


4.21(13) 


4.38(18) 


0.173(98) 




0.0630 


0.87080(78) 


0.8751(10) 


0.00430(47) 


4.30(18) 


4.51(25) 


0.21(12) 


Set-II 


0.1100 


0.29323(27) 


0.31229(52) 


0.01906(39) 


3.952(36) 


3.970(48) 


0.018(44) 




0.1020 


0.38526(30) 


0.39961(50) 


0.01435(34) 


3.905(37) 


3.936(49) 


0.031(40) 




0.0930 


0.49889(35) 


0.50990(57) 


0.01101(37) 


3.846(46) 


3.904(66) 


0.058(51) 




0.0840 


0.63258(40) 


0.64135(60) 


0.00877(36) 


3.757(54) 


3.829(76) 


0.071(54) 




0.0760 


0.77918(49) 


0.78647(75) 


0.00729(43) 


3.666(73) 


3.77(11) 


0.108(71) 




0.0700 


0.91609(54) 


0.92246(80) 


0.00637(44) 


3.571(84) 


3.69(12) 


0.118(73) 




0.0630 


1.11749(62) 


1.12294(87) 


0.00544(45) 


3.45(10) 


3.58(14) 


0.130(77) 



For the charm quark mass, the present lattice spac- 
ing a a is already well less than a^ l ~ u so that the 7f 
tuned for massless quark is applicable to the charmed 
hadron systems. On the other hand, bottom quark mass 
is not the case, and one need finer lattice spacing or larger 
anisotropy £. The present lattice would be also sufficient 
for the charmonium states, since the region m q < 0.2 
also covers the charm quark mass. Another striking fea- 
ture is that the observed £f's for heavy- light mesons are 
close to the tree level expectation. This implies that the 
deviation of £p from £ may be largely removed by the 
tree level tuning of jp- Such an approach would work 
for the spectroscopy of hadrons containing single bottom 
quark. This is a good alternative procedure to the mass 
dependent calibration using heavy-light meson, since the 
statistical error of £f from heavy-light mesons rapidly 



grows as heavy quark mass. 

Now we summarize the result for the Set-II. The heavy- 
heavy meson satisfies the relativity relation by defini- 
tion. However, the heavy-light meson dispersion rela- 
tion violates the relativity relation so that £f/£ deviates 
from unity toward smaller value as heavy quark mass 
increases. This implies that the mass dependent tun- 
ing cannot absorb the discrepancy between the £p's from 
the heavy-heavy and heavy-light meson systems. For the 
quark mass of m q < 0.2, the result for Set-II axe consis- 
tent with the result for Set-I. In that region the heavy- 
heavy meson can be successfully described, and therefore 
the result of calibration with heavy-heavy mesons is also 
valid. 

We now observe the inconsistency among the bind- 
ing energies of heavy-heavy, heavy-light, and light-light 
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FIG. 7: Fermionic anisotropy determined nonperturbatively from the dispersion relations of heavy-heavy and heavy-light 
mesons. The left and right panels show the results for the pseudoscalar and vector channels, respectively. The solid lines 
represent the tree level values according to Eq. (0). 



mesons discussed in Refs. 
probed by a quantity 



[ ITU The inconsistency is 



I : 



2SM hl -(5M hh + SM u ) 
2M 2M 



(16) 



where SM = M% — Mi, Mi and M 2 are the rest and ki- 
netic masses of the meson, respectively. The subscripts 
hh, hi, and 11 represent the quark contents of each me- 
son (h for heavy and / light quarks). We neglect the 
last term in the numerator, since the calibration of light 
quark mass region requires that SMu vanishes. Since the 
rest and kinetic quark masses cancel in each kind of mass, 
nonvanishing / represents the inconsistency in the bind- 
ing energy, namely the dynamical effect. The anomalous 
behavior of I in large kinetic mass region was first re- 
ported in Ref. (w) for the O(a) improved quark action 
on isotropic lattice. Ref. [151 explained that this behav- 
ior originates from the 0(Jap) 2 ) discretization effect in 
the heavy quarkonium system and estimate the size of / 
with the help of potential model analysis. It results in 
/ ~ —0.5 at aM 2 hi — 3.2, which is in good agreement 
with the result in Q. 

Figure ^ displays the results of I for Set-I and If for 
the pseudoscalar channel. The behavior of / is quite sim- 
ilar to that in Refs. Ill, fill, as is expected from that it 
originates from 0((aJp) 2 )error and therefore is not im- 
proved by the anisotropy. The behavior of / of two sets 
is very similar to each other. It means that the inconsis- 
tency cannot be eliminated by tuning anisotropy param- 
eters with either heavy-heavy or heavy-light system, and 
a kind of universal quantity to signal the failure in con- 
sistent description of heavy quarkonium systems. The I 
rapidly deviate from zero for the heavy quark mass larger 
than our lightest one. This is consistent with the above 
observation of £f- 

These results indicate that for a T m,Q > 0.2 the present 




FIG. 8: The inconsistency among the binding energies of 
heavy-heavy, heavy-light and light-light mesons. The hori- 
zontal axis is in spatial lattice units. 



action fails to describe the dynamics of heavy quarko- 
nium, and cannot be reinterpreted by the subtraction of 
quark mass effect. Thus we need to either give up to 
apply the present framework for such a quark mass re- 
gion, or to improve the action by incorporating higher 
order correction terms. The remaining part of this sec- 
tion therefore focus on the heavy-light meson spectrum. 



C. Heavy-light meson spectrum 

Now we turn our attention to the heavy-light meson, 
which is of our main interest. According to the heavy 
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quark expansion, the spin flipping interaction of heavy 
quark in the heavy-light systems is of 0(l/mg), and 
hence the hyperfine splitting of mesons is proportional to 
the meson mass at the leading order. In the heavy-light 
systems, the large mass of heavy quark is not important 
for the dynamics. Employing the Fcrmilab formulation, 
this character is taken into account and the the correct 
heavy quark expansion is in terms of kinetic mass 0. 
Therefore, the hyperfine splitting, which is measured as 
the difference of rest masses, are expected to be recipro- 
cally proportional to the kinetic meson mass. If this is 
the case, Set-I a,nd Set-II should show the similar behav- 
iors, up to the 0(a 2 ) systematic uncertainty. 

The hyperfine splitting m v — m PS is displayed in Fig- 
ure ^|. The horizontal axis is spin averaged kinetic me- 
son mass inverse in physical units defined through the 
hadronic radius tq. Since the most serious uncertainty 
in the kinetic mass is from the systematic uncertainty in 
£f, we estimate this error as 



5M 3 = 2M 2 \8fr\/t F , 



(17) 



where S^f = ft - ft , and hence it does not in- 
clude the statistical errors of £f and M^. On the other 
hand, the error associated to the hyperfine splitting is 
the statistical one. Although the data shows a linear 
dependence in the heavy quark mass region, there ap- 
pears small negative intercept. We regard this small dis- 
crepancy with the heavy quark expansion as the Lorentz 
symmetry breaking effect brought in the meson disper- 
sion relation and of 0(a 2 ). For fixed physical quark mass, 
this effect is expected to disappear linearly in a 2 toward 
the continuum limit. The results of Set-I and Set-II 
clearly show similar behavior, and therefore the above 
interpretation of Fermilab formulation works well up to 
the relativity breaking effect represented by <5£f- For 
more quantitative analysis, one must quantify the size of 
0(a 2 ) systematic errors and observe how they disappear 
toward the continuum limit. This is out of the scope of 
this paper. 



V. CONCLUSION 

In this paper, we investigated an applicability of the 
anisotropic lattice quark action in heavy quark mass re- 
gion, on the quenched lattice with a^ 1 ~ 1.6 GeV and 
renormalized anisotropy £ = 4. The anisotropy is ex- 
pected to extend the region in which the parameters in 
the action tuned for massless quark are applicable to pre- 
cision computations of heavy-light matrix elements. For 
this purpose, we measured the heavy-heavy and heavy- 
light meson masses and dispersion relations to observe 
the relativity breaking effect. The calculation was car- 
ried out for two sets of parameters, Set-I and Set-II. Set-I 
adopts the values at the massless limit, while in Set-IIthe 
bare anisotropy is tuned using the heavy-heavy mesons. 
Our main results are as follows. 
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FIG. 9: The hyperfine splitting of heavy-light meson versus 
the inverse of the heavy-light meson kinetic mass. The phys- 
ical scale is set by the hadronic radius ro- The error in I/M2 
is the systematic one estimated from the discrepancy of £f 
from pseudoscalar and vector mesons. 



(a) In the quark mass region a T niQ < 0.2, the observed 
fermionic anisotropy £f's are consistent for heavy-heavy 
and heavy- light mesons as well as for Set-I and Set-II pa- 
rameters within 2 % accuracy. This implies that present 
framework is applicable to both these systems even with- 
out tuning of anisotropy parameter. Beyond this region, 
the action fails to describe the heavy quarkonium state 
correctly, as expected from the tree level analysis of the 
quark dispersion relation. 

(b) The mass dependence of the renormalized 
anisotropy for the heavy-light mesons with 7^ at the 
massless limit (Set-I) is so small that one can exploit the 
massless tuning in the region of a r rriQ < 0.3 with less 
than 2% errors. This result is essentially important for 
our strategy, since it implies that the parameters tuned 
at the massless limit is directly available for this mass 
region, which already contains charm quark mass with 
the present lattice and can be extended to the bottom 
quark with development of computational resources in 
successive decade. 

(c) For a T mQ > 0.3, relativity breaking effect in the 
heavy-light mesons seems to grow as a function of the 
heavy quark mass. This is signaled by the discrepancy of 
£f's from the pseudoscalar and vector mesons, although 
the present statistics is not enough for a quantitative es- 
timate of this effect. In the scaling of hyperfine split- 
ting, we also found small discrepancy with the expecta- 
tion from the heavy quark expansion, which is considered 
as 0(a 2 ) systematic effect. It is important to quantify 
these effects and to observe how they disappear toward 
the continuum limit, for future high precision computa- 
tions with this approach. 

We therefore conclude that the anisotropic lattice 
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quark action satisfies the required characters (i)-(iii) for 
high precision calculations of heavy-light matrix elements 
mentioned in Introduction. Among above results, the 
small mass dependence of anisotropy parameter is in par- 
ticular quite encouraging for further development of the 
framework in this direction. 

One of the promising strategies is to calibrate the pa- 
rameters in the action at the massless limit and use them 
for all masses. This should work perfectly for a T rriQ < 0.3 
which is already sufficient to describe the D meson sys- 
tems. If one wants to treat B meson systems, one has 
to work with heavier quark, a T m,Q > 0.3, in which case 
the mass dependent errors cannot be neglected. How- 
ever, since the mass dependence is small, it can be inter- 
preted as an 0(a 2 ) error and can be removed by taking 
the continuum limit. Alternatively, one can also apply 
the genuine Fermilab approach for the bottom quark. In 
this case, the mass dependences of the renormalization 
coefficients are the source of systematic errors. Never- 
theless, as long as one obtains such coefficients nonper- 
turbatively in the massless limit first, and use the one- 
loop perturbation theory only to compute the mass de- 
pendent corrections, the perturbative error can be much 
better controlled compared to the Fermilab approach on 
the isotropic lattice. The nice agreement between the ob- 



served £f and the tree level expectation suggest that this 
idea is also promising in the bottom quark mass region. 

On the other hand, present level of systematic errors as 
well as the statistical one is not sufficient for these ends. 
As a most significant example, one need to eliminate the 
0(aa) error using nonperturbative renormalization tech- 
nique in the massless limit 111]], and confirm that this 
suffices also for heavy quark mass region treated in this 
paper. 
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